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ABSTRACT 

Numerical simulations carried out over the past decade suggest that the orbits of the Global 
Navigation Satellite Systems are unstable, resulting in an apparent chaotic growth of the ec¬ 
centricity. Here we show that the irregular and haphazard character of these orbits reflects a 
similar irregularity in the orbits of many celestial bodies in our Solar System. We find that 
secular resonances, involving linear combinations of the frequencies of nodal and apsidal pre¬ 
cession and the rate of regression of lunar nodes, occur in profusion so that the phase space is 
threaded by a devious stochastic web. As in all cases in the Solar System, chaos ensues where 
resonances overlap. These results may be significant for the analysis of disposal strategies for 
the four constellations in this precarious region of space. 

Key words: celestial mechanics - chaos - methods: analytical - methods: numerical - planets 
and satellites: dynamical evolution and stability — planets and satellites: general. 


1 INTRODUCTION 


Space debris—remnants of past missions, satellite explosions, and 
collisions—is a phenomenon that has existed since the beginning of 
the space age; however, its significance for space activities, in par¬ 
ticular the increasing impact risks pos ed to space systems, has been 
realised only in the past few decades llKessler & Cour-Palaisl 19781 : 
Rossi et al Jl999l : lLiou & Johnsonl2006it . The proliferation of space 
debris has motivated deeper and more fundamental a nalysis of 
the long-term evolutio n of orbits about Earth dBreiterl [200 1 bl ltl 
ICelletti & Galesll2014li . Orbital reson ances are widespread within 
this system as a whole (Hughes 1980), but particularly so amongst 
the medium-Earth orbits (MEOs) of the navigation satellites in the 
region of semimajor axes between 4 and 5 Earth radii, and a clear 
picture of their nature is of g reat importance in assessing debris mit¬ 
igation measures l lAlessi et al.l2014l) . Indeed, the discovery that the 
recommended graveyard orbits of these satellites, located several 
hundred kilometres above the operational constellations, are poten¬ 
tially unstable has led to a new paradigm in post-mission disposal— 
one that seeks to cleverly exploit these dynamical instabilities 
and the associated eccentricity growth for re-entry and destruction 
within the Earth’s atmosphe re ifjenki n & Gickll2002l : IChao & Clickl 


120041 : R OSS il boosl ; fueled ie et al.l!201 ll) . Previous studies have al 


ready noted the connection between the origin of the long-timescale 
instabilities in the MEO region and a resonance phenomenon in¬ 
volving Earth oblateness and lunisolar perturbations, yet very little 
attention has been given to a true physical explanation of the erratic 
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behaviour. To identify the source of orbital instability, we inves¬ 
tigated the main resonant structures that organise and govern the 
long-term orbital motion of navigation satellites. This paper clari¬ 
fies the fundamental role played by secular resonances and reveals 
the significance of the perturbed motion of the Moon. 


2 ANALYSIS 

Orbital resonances are ubiquitous in the Solar System and 


are harbingers for the onset of di 

mamical instability 

and 

chaos ( Duncan et al. 19951: Lecaretal. 

2001; Morbidellil 2003: 


Colwell et al J2009I : Tsiganisl20ld : Asphaug|2014l : Lithwick & Wd 

201 4) . Our knowledge of such phenomena in celestial mechanics 
comes mainly from studies on asteroid dynamics—attempts to de¬ 
scribe the stunningly complex dynamical struct ure of the as teroid 
belt since the discovery of the Kirkwood gaps jlVIorbidel 1 il 1 20021 : 
iTsiganislhoiol ). The orbital structure of the main belt is divided 
by unstable regions, narrow gaps in the distribution of the aster- 
oidal semimajor axes at 2.5 and 2.8 astronomical units (AU) from 
the Sun, where the orbital periods (or mean motions) of the as¬ 
teroids are commensurable to that of Jupiter. The outer boundary 
is marked by the Jovian 2:1 mean-motion resonance at 3.3 AU. 
Secular resonances, which involve commensurabilities amongst the 
slow frequencies of orbital precession of an asteroid and the plan¬ 
ets, define the inner edge of the main belt near 2.1 AU and demar¬ 
cate the recognised subpopulations (asteroid families and groups). 
Secular resonances also occur embedded within the libration zones 
of the prominent Jovian mean-motion resonances, and it is pre- 
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cisely the coupling of these two mai n types of resonance phenom¬ 
ena that generates w idespread chaos (Leca r et a l,. i200ll :lMorbidelli 
1 20021 : pTsi ganisl 2oToh . The creation of the Kirkwood gaps is thereby 
associated with a chaotic growth of the asteroids’ eccentricities 
to planet-crossing orbits, which leads to the removal from these 
resonant regions by collisions or close encounters with the ter¬ 
restrial planets. The same physical mechanism that sculpted the 
asteroid belt—the overlapping of resonances—is responsible for 
the emergence of chaos in the Solar System on all astronomical 
scales, from the dynamics of the trans-Neptunian_Kuiper belt and 
scattered disk populations of small bodies jDuncan et al .11 19951) to 
the a ppare n t per fectly regulated clockwork motion of the planets 
dLithwick & WulhoA 

With the advent of artificial satellites and creation of space 
debris, similar dynamical situations, every bit as varied and rich 
as those in the a steroid belt, can occur in closer proximity to ou r 
terrestrial abode l lUpton et atll 19591: lHugheslll980t iBreiteifeoOlal) . 
The dynamical environment occupied by these artificial celestial 
bodies is subject to motions that are widely separated in frequency: 
the earthly day, the lunar month, the solar year, and various preces¬ 
sion frequencies ranging from a few years to nearly 26,000 years 
for the equinoxes. A vast and hardly surveyable profusion of per¬ 
turbations originates from the gravitational action of the Sun on the 
Earth-Moon system. Primary among these irregularities is the re¬ 
gression the Moon’s line of nodes with a period of about 18.61 
years, and a progression of the lunar apsidal line with a period 
of roughly 8.85 years. The provision of frequencies in the Earth- 
Moon-Sun system gives rise to a diverse range o f complex reso - 
nant phenomena associated with orbital motions (Hughesll 1 98(1 ). 
Indeed, past efforts to identify and classify the resonances that gen¬ 
erate the orbital instability in the MEO region have been wholly ob¬ 
scured by the abundance of frequencies in the orbital and rotational 
motions of the Eart h and of the Moon, and the great disparity of 
timescales involved iDeleflie et al.lfed 1 ll : lBordovitsvna et all2014l ; 

ICelletti & Galesl2QI4l). 


While the orbits of most artificial satellites are too low to be 
affected by mean- moti on commensurabilities, even with the Earth’s 
moon l lPichmann et alJ[2014h . there exists a possibility of a some¬ 
what more exotic resonance involving a commensurability of a sec¬ 
ular pr ecession frequency with the mean motio n of the Sun or the 
Moon JCookll 1 96l : lHougl1l 1981 1 : [B^reiterll200Fbl ). Among the more 
compelling of these semi-secular commensurabilities is the evec- 
tion resonance with the Sun, where the satellite’s apsidal preces¬ 
sion rate produced by Earth’s oblateness equals the Sun’s apparent 
mean motion, so that the line of apsides follows the Sun. This is 
the strongest of the semi-secular apsidal resonances l lUpton et al .1 
1 19591 : iBreiteilfeoOlal) . and its significance has even been stressed 
in the dynamical history of the lunar orb it after its hypo thesised 
formation from an impact-generated disk jAsphauj|2014l l. as well 
as in the dynamics of the Saturnian ring system with the larger, 
more distant moons, Titan and Iapetus, playing the part of the 
Sun |Colwell_et_al. 2009). It is therefore not surprising that some 
authors (Peleflie et alj 201 ll ; iBordovitsvna et al.fl20l4) have sug¬ 
gested a possible link between these semi-secular resonances and 
the irregular eccentricity growth in the orbits of the navigation 
satellites. While it is true that such resonances exist for a wide 
variety of eccentricities and inclinations, they are only of trivial 
importance in these reg ions, a s they mainly occur at much lower 
semimajor axes tereitenhoOlal) . 

Satellite orbits resonant with respect to the perturbing influ¬ 
ence of the Earth’s gravitational fie ld are of partic u lar i mpor tance 
for the navigation satellite systems jElv & Howelll 1 99l : lEl\l20()1 : 


lRossi|[200l : ICelletti & Galesll20l4) . and such resonance problems 
have stood in the foremost rank of astrodynamical research work 
fereiteill2001al ). Undoubtedly, the most celebrated resonance phe¬ 
nomenon in artificial satellite theory is the critical inclination prob¬ 
lem, which involves a commensurability between the two degrees 
of freedom (apsidal and nodal motions) of the secular system 
1 Houghl 1 198 ll : iDelhaise & Morbidellil 1 19931 : iTremaine & Yavet j 
120141) . The extreme proximity to these critical inclinations, espe¬ 
cially in the c ase of the European Galileo an d Russian GLONASS 
constellations J R o s s il 2t )0 8l : I AI e ss i et alfeObfl) . induced researchers, 
very early, to associate the origin of the instabilities observed 
in numerical survey s to these resonances djenkin & Gickl 120021 : 
Ichao & Gickl 120041) . However, t hese inclination-dependent- only 
resonances are generally isolated Icookll 19621 : iHughesll 1980h and 
thereby exhibit a well-behaved, pendulum-like motion, which, 
when considered alone, cannot explain the noted chaotic be¬ 
haviours. The second class of geopotential resonances arises from 
the slight longitude-dependence of the geopotential, where the 
satellite’s mean motion is commensurable with the rotation of the 
Earth. Under such conditions, the longitudinal forces due to the 
tesseral harmonics continually perturb the orbit in the same sense 
and produce long-term changes in the semimajor axis and the mean 
motion. For the inclined, slightly eccentric orbits of the naviga¬ 
tion satellites, where these commensurabilities abound, the tesseral 
harmonics exhibit a multiplet structure, akin to mean-motion reso¬ 
nances, whereby the interaction of resonant harmonics can overlap 
and produce cha otic motions of the semimajor axis dElv & Howelll 
119971 : lElvl l2002h . These effects, however, are localised to a nar¬ 
row range of semimajor axis (tens of kilom etres) and are of 
much shorter period th an secular precession dElv & Howell|[l997l : 
ICelletti & Galesll2014l) : consequently, tesseral resonances will not 
significantly affect the long-term orbit evolution over timespans of 
interest. 


3 RESULTS AND DISCUSSION 


These considerations have led us to investigate the role of lu¬ 
nar secular resonances in producing chaos and instability among 
the navigation satellites. We treat only the secular resonances of 
the lunar origin, since, de spite recent implications to the contrary 
dBordovitsvna et alj2014lh the solar counterparts will require much 
longer timespans than what inte rests us he re for their destabilising 
effects to manifest themselves felvlhooj) — a consequence of the 
disparity of several orders of magnitude between the orbit preces¬ 
sion periods of the Earth and the Moon. As a basis for our cal- 
culation s we have used a convenient trigonometric series develop¬ 
ment bv lHughel d 1980h , corresponding largely to a harmonic anal¬ 
ysis of the perturbations. The lunar disturbing function is devel¬ 
oped as a Fourier series of complicated structure, whose arguments 
are combinations of the orbital phase and orientation angles of the 
satellite and the Moon, and whose coefficients depend on the size 
and shape (semimajor axes and eccentricities) of their orbits and 
the inclinations. Two considerable simplifications are possible, re¬ 
ducing this rather formidable expression in a marked degree. For 
satellites whose semimajor axis does not exceed one tenth of the 
Moon’s distance from the Earth, we can truncate the series to sec¬ 
ond order in the ratio of semimajor axes, so that the lunar poten- 
tial is approximated with sufficient accuracy by a quadrupole field 


api __ _ 

dMusenl 196ll:IDelhaise & Morbidellill 19931:iRosengren & Scheeresl 


120131 : Tremaine & Yavetzll20l4l) . To study the secular interactions, 

the short periodic terms of the disturbing function, depending on 
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the mean anomalies o f both the satellite and the Moon, can be aver¬ 
aged out l iMorbidelli[l2002h . For lack of any mean-motion or semi¬ 
secular resonances, such averaging effectively involves discarding 
all terms which depend on these fast orbital phases. Accordingly, 
the lunar disturbing function reduces to the form: 


n = 


H M a 


<(! t) 3/2 


££(-irjf„ 


(2-a)! 

(2 + m)\ 


m =0 s=0 

2 

^ ' F2,m,p(t)F2,3,l(tM)G , 2,p,2p-2(e) 

[(- 1 ) 2_S[/ 2 n ’- s (e)cos$p im!S + t/™’ s (e)cos$“ m , s ] , 
with harmonic angles 


(1) 


= (2 - 2p)uj + mQ ± sn M ± ^(s + m), (2) 

where m, s, and p are integers and the quantity K m is such that 
Kq = 1 and K m = 2 for m > 0. The semimajor axis a, eccentric¬ 
ity e, inclination i, longitude of ascending node fi, and argument of 
pericentre uj are the satellite’s orbital elements relative to the celes¬ 
tial equator. The ecliptic orbital elements of the Moon are denoted 
by subscript M, its gravitational mass by /im, and e is the Earth’s 
obliquity. The functions F and G are the modified Allan inclination 
function and Hansen coefficient, respectively; and the function U 
acc ounts for the fa ct that the Moon’s orbit is referred to the ecliptic 
(see lHughesll 198(1 for more details): 


f/ 2 m ’ ±s 


(- 1 ) 


m=ps 


(2 ± s)! 
d 2± 


(cos e/2) m±s (sin e/2) d 


AZ n±s 


{Z 2 - m (Z-l) n+ mj , 


(3) 


where Z = cos 2 e/2. Such a choice of reference planes permits us 
to consider Im as roughly constant and the variation in SIm, largely 
caused by the disturbing actio n of the Sun, as approximately linear 
with period 18.61 years jHugheslI 198(1) . The rather remarkable fact 
that the Moon’s perigee does not explicitly appear in the quadrupo- 
lar expansion of the secular problem has been k nown for some time 
l Museni n~96 lh . and belies the recent assertion dBordovitsvna et all 
l2014h that the position of the lunar perigee gives rise to resonance 
effects on navigation satellite orbits. 

A lunar secular resonance occurs when 


Ip 2 - 2 p,m ,±3 = (2 — 2p)uj + mtl ± sQ m « 0; (4) 

that is, when a specific linear combination of the secular precession 
frequencies vanishes. For the Moon, the rate of change of Qm is ap¬ 
proximately —0.053 deg/d. The oblateness precession completely 
overshadows the lunisolar effects and other perturbing influences 
so that uj and Q have an essentially linear time dependence^ with 
secular rates given by the classical expressions JCoo ll962j): 


1 This approximation is satisfactory when compared with the rigorous the¬ 
ory. Neglect of the lunisolar perturbations on the the frequencies of nodal 
and apsidal precession sets an upper limit to the radius of the orbit for which 
the theory is valid; on the other hand, the period of precession must be ap¬ 
preciably longer than a year for the double averaging procedure to be justi¬ 
fiable, which sets a lower limit to the orbital radius. For these reasons, the 
analysis is most useful in the region of semimajor axes between 4 and 5 
Earth radii. 



5 cos 2 i — 1 
(1 — e 2 ) 2 


cos i 

(1 — e 2 ) 2 ’ 


(5) 


where J 2 is the second zonal harmonic coefficient of the geopo¬ 
tential, R is the mean equatorial radius of the Earth, and n is the 
satellite’s mean motion. Equations 0 and 0 define surfaces of 
secular resonances in the space of the orbital elements a, e, and 
i. As the semimajor axis is secularly invariable after averaging, the 
commensurability condition defines curves of secular resonances in 
eccentricity and inclination phase space. The interaction of various 
lunar harmonics produces an exceedingly complicated network of 
resonances—courtesy of the Earth’s oblateness—becoming partic¬ 
ularly dense near the inclinations of the navigation satellite orbits 
(Fig.0. Each of the critical inclinations of the geopotential reso¬ 
nances (corresponding to s = (@) split into a multiplet of resonant 
curves, emanating from unity eccentricity^ Each of the resonance 
curves in Fig. |T| has its own pendulum-like structure, with a char¬ 
acteristic domain of libration. When the separation between nearby 
resonances becomes similar to their libration widths, the resonant 
critical angles of the trajectories in these overlapping regions may 
switch irregularly betw een libration and circulation; this alterna tion 
is a hallmark of chaos J Morbidel 1 ill2002l ; fLi th wick & Wj2QI4l) . 

To gain an understanding of the physical significances of these 
resonances, it it important to develop analytical models that ac¬ 
curately reflect the true nature of the resonant interactions. The 
simplest physical model for the study of the secular and resonant 
motion in MEO consists of the quadrupole order of approximation 
for the Earth, lunar, and solar gravitational potentials, but with the 
Moon’s motion represented by a precessing and rotating elliptical 
orbit. The singly-averaged model forms a nonautonomous Hamil¬ 
tonian two degrees-of-freedom system, depen ding perio dically on 
ti me through the Moons perturbed motion l IRosengren & Scheeresl 
l2013ll . Numerical simulations of this model show that the regions 
where two or more neighbouring lunar secular resonances inter¬ 
act exhibit chaos and large-scale excursions in eccentricity (Fig.0, 
in agreement wit h theo reti cal pr edictions based on the resonance 
overlap criterion fchirikovlll979ll . Chaotic diffusion is mediated by 
the web-like structure of secular resonances, which permeate the 
phase space and allows an initially circular orbit to become highly 
eccentric, as revealed in a Poincare map of the nonautonomous, 
periodically-perturbed two degrees-of-freedom system, obtained 
by sampling the flow stroboscopically once per lunar nodal pe¬ 
riod. Once a dynamical instability sets in, the subsequent evolu¬ 
tion is highly chaotic and unpredictable in detail. The resonances 
are the preferential routes for chaotic diffusion, as the trajectory 
jumps from one resonance domain to the other. In any case, the 
perigee will drive toward the ground, and the satellite is doomed 
to destruction. The chaotic zones defined by the regions of overlap¬ 
ping resonances do not preclude the existence of regular trajectories 
embedded within it. Indeed, the character of the motion depends 
sensitively upon the initial orientation angles of the satellite and 
the initial lunar node (Fig. [2]>. The second orbit (bottom) suffers 


2 iHughed j 198Ct) classifies these commensurabilities as inclination- 
dependent-only lunisolar secular resonances, under the assumption that the 
J 2 harmonic in the geopotential produces the dominant cha nge in the satel¬ 
lite’s perigee and node (see also Cook 1962; Breiter 2001a). 

3 To depict such curves up to e = 1 is only theoretical, as the satellites 
will re-enter Earth’s atmosphere when e > 1 — R/a. Nevertheless, we are 
mainly interested here in the early stages in the development of chaos. 
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Figure 1. The resonant structure of the medium-Earth orbits of the four 
navigation constellations: GLONASS (a = 4.001?), GPS (a = 4.16.R), 
BeiDou (a = 4.387?), and Galileo (a = 4.647?). Each curve repre¬ 
sents the location of the exact resonance of the form 'ip 2 — 2 p,m,±s = 
(2 — 2 p)uj + mCl =b s&m = 0, where the effects of perturbations on 
co and O other than the J 2 harmonic have been neglected. For each group 
of curves, the vertical line corresponds to the inclination-dependent-only 
lunisolar secular reson ances (s = 0), while the other curves are obtained 
fors = -2. -1.1. 2. IE1 v & Howel jl997h were the first to identify this 
devious network of resonances for the GPS constellation; yet the exact na¬ 
ture of its consequences has not been hitherto explored. These resonances 
generally form the skeletal structure of the phase space, tracing the topolog- 
ical or ganis ation o f the manifolds on which the chaotic motions take place 
iMorbidellil 120021 . Chaos i s pro duced by the overlapping of the closely- 
spaced resonant harmonics I Chirikov 1979). 

only small-amplitude variations in eccentricity and inclination, and 
is apparently regular, at least on the timescale of 500 years. Fig¬ 
ure [3] shows that the absence of resonance overlapping generally 
guarantees the local confinement of the motion. A large amount 
of simulations, in which the initial perigee and nodes were varied, 
were performed and further confirm this result: satellites with initial 
conditions in the non-overlapping region exhibit regular, bounded 
motions. For the full system one might expect that a coupling of 
the tesseral and lunar resonance phenomena would produce a slow 
diffusion in phase space, by which orbits can explore large regions 
jumping from one resonance to another. A dis cussion of t his aspect 
of chaotic diffusion has been given by others dElv & Howellll 19971 : 
lEivhood) . and only occurs over millennia timescales. 

It should be reiterated that the resonant curves presented in 
Fig. [I] show the regions in the inclination-eccentricity phase space 
for which chaotic orbits can be found and explain why chaotic or¬ 
bits manifest in these regions only. It gives, however, no informa¬ 
tion about which initial angles (to, £1, and Qm) will lead to chaos. 
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Figure 2. Dynamical evolution of two Galileo-like orbits superimposed on 
the background of complex resonant curves. The initial epoch of the sim¬ 
ulations is 26 February 1998, which determines the initial dynamical con¬ 
figuration of the Earth-Moon-Sun system, and the orbits were propagated 
for 500 years (copper-tinge sequence of discrete points) using an accurate, 
non-singular, singly-averaged model of oblateness and lunisolar perturba¬ 
tions, develop previously by one of the authors ^Rosengren & Scheeresl 
hoiah . Both representative Galileo orbits initially have a = 29, 600 km, 
e = 0.001, i = 56°, u> = 30°, but their orbital planes are separated by 
120° with f i = 240° (top) and S2 = 120° (bottom), respectively. We con¬ 
sider the discrete dynamics of the system by taking a subsequent snapshot 
of the motion every lunar nodal period (a stroboscopic map), which gives 
a faithful picture of the general character of the trajectories. The motions 
with initial conditions in the overlapping regions can be both chaotic and 
regular, depending on the initial orientation angles. 

The structure of the whole phase space is far from being fully un¬ 
derstood and clearly warrants greater theoretical development. Fu¬ 
ture work will present more complicated maps, which allow for 
the global visualisation of the geometrical organisation and coexis¬ 
tence of chaotic and regular motion. 

4 CONCLUSION 

In the inner Solar System, overlapping secular resonances have 
been identifi ed as the origin o f chaos in the orb i ts of the terres- 
trial planets iLecar et all 1200ll : iMorbidellil 120021 : iLithwick & Wul 
120141) . We have demonstrated that the same underlying dynamical 
mechanism responsible for the eventual destabilisation of Mercury, 
and recently proposed to explain the orbital architecture of extra¬ 
solar planetary systems dLithwick & Wul2Q14l) . is at the heart of the 
orbital instabilities of seemingly more mundane celestial bodies— 
the Earth's navigation satellites. The occurrence and nature of the 
secular resonances driving these dynamics depend chiefly on one 
aspect of the Moon’s perturbed motion, namely, the regression of 
the line of nodes. The decisive significance of this fact on space de¬ 
bris mitigation will be emphasised in a later paper an d its rel ation 
to the design of disposal strategies in MEO dAlessi et alj2014l) will 
be formulated there more precisely. The precarious state of these 
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Figure 3. Dynamical evolution of five orbits with initial inclinations and 
eccentricities in non-overlapping resonance regions superimposed on the 
background of secular resonant curves. The initial epoch of the simu¬ 
lations is 26 February 1998, and the orbits were propagated for 500 
years (copper-tinge sequence of discrete points) using an accurate, non- 
singular, singly-averaged mod el of oblateness and lunisolar perturbations 
jRosengren & Scheered l2013h . All orbits initially have a = 29, 600 km, 
uj = 30°, but their orbital planes are separated by 120° with = 240° 
(top) and = 120° (bottom), respectively. We consider the discrete dy¬ 
namics of the system by taking a subsequent snapshot of the motion every 
lunar nodal period (a stroboscopic map). The motions with initial conditions 
in the non-overlapping regions are generally regular; chaotic motion can 
only exist close to the separatrices of the isolated resonances iMorbidelhl 

12002I) . 


constellations, perched on the threshold of instability, makes it un¬ 
derstandable why all efforts to define stable graveyard orbits, espe¬ 
cially in the case of Galileo, were bound to fail; the region is far too 
complex to allow of an adoption of the simple geosynchronous dis¬ 
posal strategy. A full understanding of the nature and consequences 
of the chaos in these environments would have certainly helped in 
the early design phases of the constellations. 

Active debris removal, apart from the daunting obsta- 
cles in engineering that this represents, is currently seen 
dLiou & Johnsonl 1200611 as the only viable option to prevent the 
self-gen erating Kessler syndrome phe nomenon (collisional cascad¬ 
ing; see lKessler & Cour-Palaisiri978h from occurring in low-Earth 
orbit, the most densely populated orbital environment dRossi et alj 
1 1999h . But as we are still remarkably ignorant of the locations and 
consequences of most resonances in near-Earth space, such drastic 
measures may require a reassessment. This concerns particularly 
the question as to whether strong instabilities exist, whose desta¬ 
bilising effects occur on decadal timescales, that can be exploited 
to effectively clear these regions of space from any future collision 
hazard. Indeed, the process of dynamical clearing of resonant or¬ 
bits is well illustrated by the paucity of asteroids observed in the 
Kirkwood gaps dTsieanisll2010h . 
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